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Abstract 

Based on a sequence of numerical computations, a conjecture is presented re- 
garding the class of functions H{x;a) = exp(a) — (1 + a/x)^, and the open 
problem of determining the values of a for which the functions are completely 
monotonic with respect to x. The critical value of a is determined here to 
sufficient accuracy to show that it is not a simple symbolic quantity. 



1. Introduction 

A completely monotonic (CM) function is an infinitely differentiable function 
/ : (0, oo) M whose derivatives satisfy 

Vx > 0,Vn e N, (-l)"/(")(a;) > . 

The definition is due to HausdorfF and a collection of important proper- 
ties can be found in 0. Alzer and Berg f]\ state the following open problem. 
Consider the function 

H{x;a) =e'' ~ (l + ^y . (1) 

For what values of the parameter a is H{x; a) completely monotonic (CM)? 

A problem that is superficially similar to this problem has a remarkably 
simple solution jjy . The function 

J{x;a,b)= + -e" (2) 

is CM if and only if a < 26. Although H{x;a) — —J{x; a,0), the problems 
for J and H are essentially different, because a completely monotonic function 
must be positive and decreasing. A more completely descriptive name might 
be 'completely monotonically decreasing function', but the shorter name is now 
standard. In spite of the beautiful result for J, a consequence of the conjecture 
advanced here is that no similarly simple result applies to H{x] a). 

Alzer and Berg [l| show that H{x; 1) is CM. Moreover, Berg Q has shown, 
through an equivalent problem, that H{x; 3) is not CM. Therefore, it has been 
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Figure 1: (a) f(x, a, 2) for three values of a. The function f{x, 2.9, 2) is monotonic, f{x, 3.4, 2) 
is not monotonic, and /(x, 3.138, 2) is the transitional case for n = 2. (b) The right plot 
shows for the critical value a = 3.138 the successive derivatives /(a;, 3.138, 2), /(x, 3.138, 3) 
and /(x, 3.138, 4). 

conjectured that there exists a value Oc such that H{x; a) is completely mono- 
tonic for all a < flc and not for a > Oc. We obtain here experimentally an 
estimate for its value, namely Ac « 2.299656443. 

In the field of experimental mathematics, there are several tools for identify- 
ing a real number with an expression containing known mathematical constants. 
We have applied these tools to Oc, but without success. 



2. Heuristic description of method 

In this section we give a graphical and heuristic description of the method 
we use. For convenience, we define the function 



(3) 



and consider its properties for various fixed n. 

We start with the case n = 2. Figure [TJa) shows plots of f{x,a,2), or 
equivalently the second derivative of H. For a = 2.9, the function is clearly 
monotonic decreasing; for a = 3.4, the function is clearly not monotonic. The 
transition occurs for a = Of ~ 3.138, when /(a;, af,2) has an infiexion point at 
X = X( ^ 0.913, as shown. 

The transitional case is further considered in figure [rfb), where the second, 
third and fourth derivatives are plotted. As could be anticipated from the shape 
of /(x, af,2), its derivative f{x,a{,3) touches the x-axis at x = Xf. Further, 
/(x, af,4) crosses the axis at the same point. From these observations, it is 
obvious that the pair of values Of , Xf satisfy the equations 



f{xi,ai,3) 



(4) 
(5) 
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Now consider the case n — 3. The above series of observations can be 
repeated for f{x,a,3). Again there is a range of values of a for which the 
function is monotonic decreasing, and there is a critical pair af,X{ such that 
/(a;f,af,3) is an inflexion point. Solving the pair of equations /(a;f,af,4) = 
/(xf,at,5) = gives at w 2.903, Xf w 1.344. In general, the function f{x,a,n) 
will be monotonic decreasing for a 'small enough'. From Alzer and Berg 
a = 1 is always such a value. The change to non-monotonic behaviour will 
occur at a value af(n) when f{x,a,n) has an inflexion point at Xf(ri). The 
values of af(n) and Xi(n) can be determined from the equations 

/(xf , af , n + l) = f{x[, Of , n + 2) = . 

We next consider how af(n) depends upon n. We have seen already that 
af(2) « 3.138, and af(3) w 2.903. For all n > 2, a[{n + 1) < ai{n). This is 
so because if f{x,a{{n),n) has an inflexion point, then f(x,a{(n),n + 1) must 
be non-monotonic. Therefore, if one thinks of a as increasing from a — 1, for 
which value all functions are monotonic, then f{x,a,n + 1) must have become 
non-monotonic at a lower value of a than that for which f{x^a,n) becomes 
non-monotonic. Thus the sequence {af(n)} must be monotonically decreasing 
and it is bounded below by 1 from the results in [1]. Therefore the limit Oc = 
lim„_^oo Of ("■) exists and is the critical value of a such that H{x;a) is completely 
monotonic for a < Oc. 

3. Numerical method 

In view of the above observations, we need to compute af (n) for a selection 
of values of n and then extrapolate to compute the limit n oo. In order 
to obtain a value of Oc that is accurate enough to be submitted to tools such 
as Plouffe's inverter [HI, we need to compute derivatives up to the order of 
10^. Hence, an efficient method must be found to evaluate f{x,a,n) for large 
values of n. Symbolic differentiation becomes impossible after at most n = 50, 
depending upon the memory available on the computer used, implying the need 
for a numerical scheme. 

By introducing the notation 

'={'+x) ' 
we can write H(x; a) ~ exp(a) — h, and obtain 

/(a;,a,0) = H{x; a) ^ e"" - h , (6) 
J'/ IN dh / x + a a \ 

= d-x-r^-^a)'- 

Leibnitz's rule now allows us to compute recursively. 

f(x,a,n) = (-1)' 

= (-1)' 



.1 d" 



dx" 

n-l 
fe=0 



jf{x,a, 1) 

1\ d'^h d" 
dx'^ dx^ 



-l-k 



-i-fe 



In- 



.(8) 
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The second term in the sum can be computed explicitly by introducing the 
following function. 



I - In + -S_ ^ 




(x-\-a) 



+ 



(2;+a)" + i ' 



n = , 
n > 1. 



Thus the computational form of (jH)) becomes 



/(a;,a,n) = ^f ^ jf{x,a,k)g{. 




x,a,n — I — k) . 



(9) 



Although it is possible to compute with this formula symbolically, expression 
swell prevents this approach from being useful for anything other than for error 
checking. If we assign numerical values to x and a and compute /(a;,a, 0) 
numerically, then ([9]) can evaluate derivatives numerically to orders 10^ or more. 

Once the functions f{x,a,n) can be computed efficiently, the equations 
f{x{,a{,n + 1) = 0, f{x{,af,n + 2) = can be solved to find the transition 
point in f(x,a,n). For small values of n, any method, for example Maple's 
f solve, can be used. For large n, extended precision must be used, owing to 
loss of precision through the accumulaton of cancellation and rounding errors. 
For n greater than 10'* the precision loss amounts to about 15 decimal dig- 
its. The equations were solved using bivariate Newton iteration. The required 
derivatives with respect to a can be computed using a method analogous to that 
above. In addition to cancelation errors, there is another reason for requiring 
extended precision. Because {xf, at) is a repeated root of f{xf, af , n + 1), New- 
ton iteration converges slowly and a nearly singular matrix must be inverted. 
Therefore 60 decimal digits were used in the computations. For smaller n, Maple 
was used, but became too slow for larger n, and Bailey's ARPREC Q and Gnu 
GMP were used. 

4. Results 

Tabled] displays critical values computed between n — 1000 and n — 100000. 
Each entry gives the smallest value of a such that f{xf, a,n + 1) < 0. That is, 
for each entry, f{xf, af — 10"-^'*, n -I- 1) > 0. 

To extrapolate to n cxd, we examine Of as a function of n^^. Figure [2] 
shows a plot of af{n) against n^^. We conjecture that af obeys a relation 



By fitting a quadratic in n^^ using the data points n — 40000, 50000, 100000, 
we obtain Oc ~ 2.29965644325, where the error can be expected to be O(10~^^. 
This value was given to Plouffe's inverter [5], the inverse symbolic calculator p, 
and Maple's identify command. No symbolic quantity matched all digits, and 
the closest symbolic quantities offered no immediate inspiration for analytic 
investigations. 



flf = ac + - + ^ + 0(n-3). 



(10) 
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n 






1000 


2.30183958971854 


436.380167908055 


2000 


2.30074838075010 


872.743540008136 


4000 


2.30020250313093 


1745.47034071250 


5000 


2.30009330574014 


2181.83374672043 


10000 


2.29987488908353 


4363.65078807689 


12500 


2.29983120225165 


5454.55931101894 


16000 


2.29979297531646 


6981.83124393025 


20000 


2.29976566981568 


8727.28488210931 


40000 


2.29971105744643 


17454.5530758350 


50000 


2.29970013475374 


21818.1871732641 


100000 


2.29967828914950 


43636.3576615316 



Table 1: Table of critical values. 



5. Closing remarks 

The simple solution to problem ([2]) has led to several attempts to find an 
equally simple solution for ([T]) . The numerical evidence presented suggests that 
such attempts will continue to fail. The numerical data obtained here does, 
however, suggest a number of simple properties which might help in an analytical 
solution to the problem. The dependence of X{ (n) on n suggests a strongly linear 
correlation between X{ and n. A similar correlation between ac{n) and has 
been noted above. If these numerical observations can be explained analytically, 
then a full solution might follow. 
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Figure 2: A plot of af against n~^. 
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